First-principles study of lattice instabilities in the ferromagnetic martensite Ni2MnGa 
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The phonon dispersion relations and elastic constants for ferromagnetic Ni2MnGa in the cubic 
and tetragonally distorted Heusler structures are computed using density-functional and density- 
functional perturbation theory within the spin-polarized generalized-gradient approximation. For 
0.9 < c/a < 1.06, the TA2 tranverse acoustic branch along [110] and symmetry-related directions 
displays a dynamical instability at a wavevector that depends on c/a. Through examination of the 
Fermi-surface nesting and electron-phonon coupling, this is identified as a Kohn anomaly. In the 
parent cubic phase the computed tetragonal shear elastic constant, C'=(Cii— Ci2)/2, is close to 
zero, indicating a marginal elastic instability towards a uniform tetragonal distortion. We conclude 
that the cubic Heusler structure is unstable against a family of energy-lowering distortions produced 
by the coupling between a uniform tetragonal distortion and the corresponding [110] modulation. 
The computed relation between the c/a ratio and the modulation wavevector is in excellent agree- 
ment with structural data on the premartensitic (c/a — 1) and martensitic (c/a = 0.94) phases of 
Ni2MnGa. 
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I. INTRODUCTION 

Ferromagnetic shape-memory alloys displaying large 
magnctic-ficld-induccd strain have recently emerged as 
a new class of active materials, very promising for actu- 
ator and sensor applications. The largest known mag- 
netostrain effects have been observed in Ni2MnGa-based 
Heusler alloys, where up to 6% and 9.5% strains have 
been induced by a magnetic field less than 1 Tcsla.^'^ 

The magnetic shape-memory behavior is closely linked 
to the occurrence of a martensitic transformation in con- 
juction with a strong magnetocrystallinc anisotropy of 
the low-temperature martensitic phase. For a deeper un- 
derstanding of the magnetic shape-memory mechanism, a 
microscopic explanation for the origin of the martensitic 
transformation and the magnetocrystallinc anisotropy is 
clearly needed. Towards this goal, we have used ab ini- 
tio techniques to investigate the origin of the martensitic 
transformation. 

Ni2MnGa is the most intensively studied of the rel- 
atively few known ferromagnetic shape-memory materi- 
als. A number of thermal and stress- induced martensitic 
transformations have been observed in Ni2MnGa-bascd 
alloys, and both transition temperatures and crystallo- 
graphic structures are quite sensitive to alloy composi- 
tion. 

In the stoichiometric alloy, Ni2MnGa, which is fer- 
romagnetic below Tc « 380 K, two thermally- induced 
phase transitions have been observed. From the high- 
temperature cubic Heusler structure, a premartensitic 
phase transformation to a modulated cubic structure oc- 
curs below TpM ~ 260 K,^''' followed by a martensitic 
transformation to a modulated tetragonal structure be- 
low Tm ~ 220 K.^ The low temperature martensitic 
phase has a tetragonal structure (c/a=0.94, a=5.90A) 
with a superimposed incommensurate modulation along 
the [110] direction, consisting of a shuffling of (110) 
planes in the [TlO] direction with a periodicity of al- 
most 5 interplanar distances, which corresponds to a 
wave vector qM « ^{OAS, 0.43, 0) (Rcfs. 4, 6). The 
low-tcmpcrature martensitic phase has a high magne- 
tocrystallinc anisotropy,'' making it useful for applica- 
tions. At a phenomenological level, the coupling between 
strain, modulation, and magnetization has recently been 
described in a Landau- theory framework.^ 

Some clues to the microscopic origin of the martensitic 
transitions are provided by the softening of particular 
phonons and elastic constants. Inelastic neutron scat- 
tering experiments on the high-temperature phase found 
a significant, though incomplete, softening in the TA2 
phonon branch along the [110] direction, at a wave vec- 
tor qpM « ^(0.33, 0.33, 0) (Ref. 3). The phonon soft- 
ening has been found to correspond to a premartensitic 
phase occurring between the high-temperature cubic and 
low-temperature martensitic structures. The premarten- 
sitic phase has a cubic structure with a superimposed 
[110]-transverse shuffling modulation, analogous to the 



modulation of the martensitic phase but with a different 
periodicity of almost 6 interplanar distances.'^ While the 
PM phase is anticipated by a precursor phonon softening 
at qpM, no phonon softening at qm has been observed 
above Tm. 

The aim of this study is to provide a unified ex- 
planation for the microscopic origin of the rich vari- 
ety of phase transitions and modulated structures oc- 
curring in this ferromagnetic shape-memory alloy. To 
this end we have performed a first-principles study of the 
phonon dispersions and lattice instabilities of ferromag- 
netic Ni2MnGa and of their dependence upon uniform 
tetragonal strain. In our calculations we use the spin- 
polarized generalized-gradient approximation (cr-GGA) 
recently developed within the density-functional pertur- 
bation theory (DFPT) formalism.® This approach allows 
us to obtain an accurate description of the structural, 
magnetic, and vibrational properties of Ni2MnGa. 

The paper is organized as follows. The details of the 
computational method are given in section II. In section 
III we present the results for the cubic Heusler structure. 
The crystal structure and magnetization are discussed in 
Section III. A. Section III.B is devoted to the phonon dis- 
persion, with particular attention to the phonon anomaly 
and related dynamical instability. The origin of the 
phonon anomaly is discussed in Section III.C. The com- 
puted elastic constants are presented and discussed in 
Section III.D. In Section IV we investigate the depen- 
dence of the phonon dispersions and lattice instabilities 
on uniform tetragonal strain. In Section V the main re- 
sults are summarized and the microscopic origin of the 
observed phase transitions explained. 

II. COMPUTATIONAL METHOD 

Our calculations have been performed within the 
framework of density-functional and density-functional- 
perturbation theory (DFT and DFPT). In particular, 
the vibrational properties have been computed using a 
recent implementation of ultrasoft pseudopotentials into 
DFPT.^° We used ultrasoft pseudopotentials^^ for Ni and 
Mn, freezing the 3s and 3p core electrons and treating the 
3d and 4s states as valence levels. ^^'^^ The nonlinear core 
correction is used to account for the overlap between the 
core and the valence charges.^** For the Ga atom we used 
a norm- conserving pseudopotential, which also includes 
a nonlinear core correction, treating the 4s and Ap states 
as valence levels. To describe the effects of exchange 
and correlation, we used the Perdew-Burke-Ernzerhof^^ 
functional with the spin-polarized generalized gradient 
correction (cr-GGA), recently implemented within the 
DFPT formalism in the PWscf codc.^ For comparison, we 
have performed calculations within the local spin density 
approximation (LSDA) using the Perdew-Zunger param- 
eterization of the exchange and correlation energy. The 
plane-wave basis set had a kinetic-energy cutoff of 25 
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Ry. The augmentation charges, required by the use of 
ultrasoft pseudopotentials, were expanded with an en- 
ergy cutoff of 450 Ry. The Brillouin zone integration was 
performed using the smearing technique, suitable for 
metaUic systems. The structural properties and most 
of the phonon frequencies are well converged using a 
first-order smearing function with a smearing parameter 
cr=0.03 Ry and an fee (6, 6, 6) Monkhorst-Pack grid,^'' 
yielding 28 k-points in the irreducible wedge of the Bril- 
louin zone (IBZ) for the cubic structure. To obtain an ac- 
curate description of the anomalous TA2 phonon branch, 
a smaller smearing parameter is necessary and, conse- 
quently, a finer k-point sampling. An accuracy to within 
a few cm^^ has been obtained using cr=0.01 Ry and an 
fee (10, 10, 10) k-point grid, yielding 110 k-points in the 
IBZ, for the cubic structure. For the tetragonal struc- 
ture, an (8, 8, 8) k-point grid yielding 144 k-points in 
the IBZ of the face-centered orthorhombic unit cell was 
used. 

To compute the full phonon dispersions of the cu- 
bic structure, we computed the interatomic force con- 
stants by Fourier transformation of the dynamical matri- 
ces computed on a (6, 6, 6) q-point grid in the fee BZ. 
The phonon dispersions along the [110] ([Oil]) direction 
were obtained by interpolating the dynamical matrices 
computed on a finer mesh of 24 q-points between F and 
the shortest reciprocal-space vector parallel to the [110] 
([Oil]) direction. An even denser mesh of q-points has 
been used in the proximity of the anomaly for a more 
accurate interpolation of the TA2 branch. 

III. CUBIC NiaMnGa 

In this section, we report the results of calculations for 
Ni2MnGa in the fee L2i Heusler structure (see figure 1). 

A. Crystal structure and magnetization 

In table I are given the minimum-energy lattice param- 
eter, ao, the bulk modulus, Bq, and the magnetic moment 
per unit cell, /ig, computed within the spin-polarized cr- 
GGA. Note that to obtain convergence to within O.Ol/is, 
Ho has been computed on a (10, 10, 10) k-point grid 
with £7=0.01 Ry. These results agree well with those 
of all-electron calculations (FLAPW) performed using 
the cr-GGA."'^^ The theoretical results arc also in very 
good agreement with the experimental data; the theoret- 
ical lattice parameter obtained with cr-GGA is equal, to 
within theoretical precision, to the experimental value. 
For comparison, the LSDA results are included in ta- 
ble I. The LSDA lattice parameter is 2.5% smaller than 
the experimental value. This underestimate of the lat- 
tice parameter is correlated with an increase in the bulk 
modulus to 38% higher than the experimental value. 



B. Phonon dispersion relation 

The phonon dispersions of the ferromagnetic cubic 
structure (oq = 11.03 a.u.) have been computed, using 
cr-GGA, for q along high-symmetry lines in the first Bril- 
louin zone, as shown in figure 2. The solid lines indicate 
the computed phonon dispersion converged to within a 
few wave numbers. The theoretical phonon dispersion 
curves are in excellent agreement with the available in- 
elastic neutron scattering data.'*'^^ 

To evaluate the accuracy of the theoretical method, the 
phonon frequencies computed using different approxima- 
tions for the exchange and correlation energy, at the zone 
boundary X point, are given in table II. Since the experi- 
mental frequency of the anomalous (CC0)-aco\istic branch 
is not strongly temperature dependent at this q-point, it 
is a good reference for evaluating the accuracy of the 
theoretical method. The experimental frequency for the 
transverse acoustic mode X5 is in very good agreement 
with the phonon frequency computed using the cr-GGA, 
with the computed value being about 6% softer. For com- 
parison, we have also done the computation with LSDA. 
Calculations for the cubic structure with the LSDA lat- 
tice constant yield frequencies about 10% harder than 
those computed using the cr-GGA, except for the lowest 
X5 mode which is 7% softer, increasing the discrepancy 
with experiment. In addition, the ordering of the modes 
in LSDA is different. Specifically, the X4 mode is softer 
than the nearby X5 and X5 modes. We have also com- 
puted the phonon frequencies within the LSDA but fixing 
the lattice parameter to the equilibrium value computed 
in the cr-GGA. The result is that all the modes soften, 
so that the frequencies are about 10% smaller than with 
cr-GGA. This softening results in an even greater discrep- 
ancy with the experimental Xg mode (20%). We conclude 
that the better accuracy of cr-GGA is not merely an effect 
of the more accurate value of the equilibrium lattice pa- 
rameter. All calculations reported below were performed 
with fj-GGA. 

The most striking feature of the phonon dispersion 
relation in figure 2 is the anomalous dip in the low- 
est branch of the transverse acoustic modes (TA2) along 
[110]. The minimum at the incommensurate wave vec- 
tor qo=^(Co, Co, 0), with Co = 0.34, is at imaginary 
frequency. Thus, the crystal is dynamically unstable to 
the lattice distortion corresponding to the eigenvector of 
this mode, which consists of a nearly rigid displacement 
of the (110) atomic planes along the [110] direction with 
a modulation period of slightly less than 6 interplanar 
distances along the [110] direction. This energy- lowering 
distortion can be specified by an amplitude, u, and a 
phase, (/), 

Um = u cos(rnCo7r -|- 0), (1) 

where u„i is the displacement of the m*'' (110) atomic 
plane along the [110] direction. Due to the cubic symme- 
try, there are 12 equivalent anomalies along the <110> 
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directions which give rise to six different but crystallo- 
grafically equivalent lattice modulations u^. 

The phonon anomaly is shown in more detail in fig- 
ure 3, where the theoretical dispersion is compared with 
inelastic neutron scattering data taken at two different 
temperatures: T=370 K (squares) and T=250 K (trian- 
gles).^ The experimental data display an anomaly in the 
TA2 branch that corresponds to the anomaly predicted 
by our theoretical dispersion. The wave vector of the 
experimentally-observed anomaly is Cpm ~ 0.33,'^ in ex- 
cellent agreement with the theoretical value. For those 
modes that do not have a strong dependence upon tem- 
perature, such as the LA branch and the TA2 modes away 
from the anomaly > 0.45), the agreement between 
theory and experiment is very good. To compare theory 
and experiment for modes that are strongly temperature 
dependent, we need to extrapolate the experimental val- 
ues to T = K. In accordance with the soft mode the- 
ory, {hu})"^ is experimentally observed to decrease linearly 
with temperature above TpM=260 K. Extrapolation to 
T = K gives an imaginary frequency of 30i cm~^, which 
is in reasonable agreement with the computed value. 

Examination of the phonon dispersion throughout the 
entire BZ shows that imaginary frequencies leading to dy- 
namical instabilities occur only in a very localized region 
in q-space. This is in agreement with the experimen- 
tal dispersion measured in the direction perpendicular 
to qo, q=qo+(Ci ~Ci O).** The contour plot in figure 4 
shows that the lattice instability is confined to a "drop- 
shaped" small region in reciprocal space, located along 
the [110] direction and centered at the critical wave vec- 
tor qo, where the dominant instability occurs. The lo- 
calized nature of the anomaly in q-space is a signature of 
its electronic origin, as will be discussed in detail in the 
next section. 



C. Origin of the anomaly 

Screening due to electron-phonon coupling involving 
electronic states near the Fermi level in metals can give 
rise to anomalous dips in the phonon dispersion, called 
Kohn anomalies. The occurrence of these anomalies de- 
pends mainly on the geometry of the Fermi surface, as 
well as on the q-dependence of electron-phonon matrix 
elements. If the Fermi surface has flat portions with nest- 
ing vector qo, there generally will be very strong screen- 
ing of the potential perturbation due to atomic displace- 
ments at that wavevector, leading to a pronounced soft- 
ening highly localized in q-space. The effect of Fermi 
surface geometry can be quantified by calculation of the 
generalized susceptibility, as in Rcf. 20. 

In cubic Ni2MnGa, the Fermi level crosses both minor- 
ity and majority spin bands. Of the two Fermi surfaces, 
only the one for the minority spin bands, plotted in fig- 
ure 5, shows obvious nesting features. In particular, a 
large fraction of the opposite sides of the "pipes," running 



along the faces of the cube and crossing at their center, 
are nested by wavevectors of the form ((q, k,0)^ with Co 
— 0.34 and < fc < 0.34 and the cubic-symmetry-related 
equivalent wavectors. Along [110], this corresponds to 
the critical wave vector of the anomaly. It should be 
noted that this differs from the position of the nesting 
wavevector and peak in generalized susceptibility along 
[110] found in Refs. 20 and 21. This difference might be 
attributable to the differences in method, in particular, 
the difference in choice of density functional. [110] ap- 
pears to be a direction for which the anomaly is strongest. 
For example, although there is nesting along [100] as well, 
there is no sign of an anomaly in the phonon dispersion. 
This can be attributed to the vanishing of the relevant 
electron-phonon matrix elements. 

The strength of a Kohn anomaly is expected to be very 
sensitive to electronic temperature, with an increasing 
temperature reducing the sharpness of the Fermi surface 
and thus weakening the anomaly. In our calculations, 
the sensitivity to electronic temperature can be directly 
investigated, as the smearing parameter a plays the role 
of a fictitious electronic temperature. In table III we 
show the dependence on a of those modes at the critical 
wave vector, qo, that have the same symmetry as the soft 
TA2 branch (i.e. the three modes of S3 symmetry). The 
anomalous TA2 mode, £3(1), is much more sensitive to 
the fictitious electronic temperature than the other two 
modes. All the frequencies arc well converged for (7=0.01 
Ry. The dependence of the phonon anomaly upon the 
electronic temperature is also shown in figure 2, where 
the dispersion of the anomalous TA2 branch computed 
with cr=0.03 Ry (dashed line) is compared with the fully 
converged calculation for fT=0.01 Ry (solid line). Only 
the modes with wavevectors close to the anomaly, (Co — 
0.16) < C < (Co + 0.16), are affected by the change in 
the electronic temperature. The anomaly is smoothed 
out as the fictitious electronic temperature is increased. 
Thus, we conclude that the computed dip is indeed a 
Kohn anomaly. 

D. Elastic constants 

The computed elastic constants for the ferromagnetic 
cubic structure arc reported in table IV. In a first- 
principles framework, the standard approach for obtain- 
ing elastic constants is to compute the stress tensor for 
a selected set of small strains. To obtain convergence 
at the level of a few GPa, we have used very dense k- 
point meshes (up to (14 14 14), corresponding to 10976 
points in the BZ). To compute Cn and C12 the cubic 
crystal is distorted by the tetragonal deformation: £22=5, 
^xx=£yy=Q- For small deformations the stress, (Jy , is lin- 
ear with e: axx=(^vy=Ci2e, and azz=Ciie. The values so 
obtained for the elastic constants Cn, C12, and the shear 
modulus C'=(Cii— Ci2)/2 are indicated as "theory 1" in 
table IV. 
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The small magnitude of C requires additional atten- 
tion. We performed an independent computation by con- 
sidering the tetragonal deformation exx=£yy=e, €22=— 2e. 
For small e we have that axx='^yy='^C'e and 0-^2=— 4C'e. 
This direct calculation gives C'=(2±2) GPa, in good 
agreement with the value obtained as a difference of Cn 
and C12 in theory 1. The near-zero value of C implies 
that there is almost no energy cost for a small tetragonal 
shear distortion of the type 2exx^'2'^yy=^^zz- Evidence 
of this marginal elastic instability towards small volume- 
preserving tetragonal distortions was also found in Ref. 
22, where it was shown that the energy surface as a func- 
tion of c/a is remarkably flat; specifically, the change in 
energy associated with varying the c/a ratio in the range 
of values between 0.97 and 1.01 is almost zero to within 
numerical accuracy. 

In table IV (theory 2) we also show the values of the 
elastic constants estimated from the slope of the long 
wavelength acoustic modes along the [110] direction. The 
elastic constants C44, C, and Cl = (Cii+Ci2+2C44)/2 
correspond to the TAi, TA2, and LA modes, respectively. 
These determine the values of Cn and C12 given in the 
table. The evaluation of the elastic constants from the 
phonon dispersion is less accurate, especially for C where 
the corresponding TA2 branch deviates from a linear be- 
havior already at very small q. 

The elastic constants in the high-temperature cu- 
bic phase have been experimentally determined between 
room temperature and the pre-martensitic structural 
phase transition at T = 260 K by measuring the veloc- 
ity of ultrasonic waves. These experiments show that 
the elastic constants Cl and Cn are almost temperature 
independent. In contrast, the transverse elastic constants 
C44 and C exhibit an anomalous behavior with cooling, 
softening as the pre-martensitic phase transition temper- 
ature is approached. While the softening is small for C44, 
it is more dramatic for C", which decreases by 60% from 
room temperature to the transition. We find good agree- 
ment between theory and room temperature measure- 
ments for Cl and Cn, which are not strongly modified 
by finite temperature. The pronounced softening of C 
is also consistent with the very small zero-temperature 
theoretical value, C ~ 0. 

IV. TETRAGONAL NijMnGa 

The identification of mistablc phonons at qo = ^(Co, 
Co, 0), with Co = 0.34, in the cubic Heusler structure 
leads naturally to an understanding of the transition to 
the premartensitic cubic-modulated phase with decreas- 
ing temperature. To explore the subsequent transition to 
the low-temperature martensitic phase with c/a = 0.94 
and modulation wavevector Co = 0.43, we extended the 
calculations of phonon dispersion to Heusler structures 
with uniform volume-preserving tetragonal strains rang- 
ing from c/a = 0.88 to c/a = 1.06, with particular atten- 



tion to the TA modes along the <110> directions. 

A. Tetragonal structure with c/a— 0.94 

We first consider in detail the volume-preserving 
tetragonal distortion with c/a=0.94, corresponding to 
the low temperature martensitic phase. The phonon dis- 
persions computed along the inequivalent [110] and [Oil] 
directions of the tetragonal BZ are shown in figure 6. 
The only mode substantially affected by the tetragonal 
distortion is the anomalous TA2 mode. We find an over- 
all softening (hardening) of the TA2 branch in the [110] 
([Oil]) direction. The anomaly is particularly affected 
by the tetragonal distortion. Along the <110> direc- 
tions, which are perpendicular to the c axis, the phonon 
anomaly is more pronounced and is shifted to a larger 
q-vector (C = 0.43) than in the cubic structure. This 
wavevector is in excellent agreement with the long pe- 
riod modulation observed experimentally in the tetrag- 
onal martensitic phase. In contrast, along the <011> 
directions, the anomaly has almost completely disap- 
peared; the entire branch is stable and there is only a 
very small wiggle in the TA2 branch. Thus, the entropy 
from low-frequency phonons should be roughly a factor 
of three less than in the premartensitic phase, which may 
explain why the premartensitic phase is more favorable 
at higher temperatures. 

As in the cubic phase, the anomaly can be associated 
with features of the Fermi surface. In figure 7 we show 
the minority-spin Fermi surface computed for Ni2MnGa 
in the tetragonal structure. It is related to the Fermi 
surface of the cubic structure with two major differences: 
(i) there arc no flat surfaces perpendicular to the (001) 
direction, and (ii) the "pipes" running along the faces 
of the tetragonal BZ are wider than the corresponding 
featiires in the cubic Fermi surface. As a consequence of 
these changes, induced by the tetragonal distortion, there 
is significant nesting only in the <110> directions per- 
pendicular to the c axis and the edges of the "pipes" are 
nested by the larger {( C, 0) critical wave- vector, C, = 0.43. 
This explains why in the tetragonal structure there are 
pronounced Kohn anomalies only in the <110> direc- 
tions, at a larger wave vector than in the cubic structure. 

B. Evolution of phonon anomaly with tetragonal 
deformation 

To understand better the lattice instabilities and re- 
lated structural energetics of Ni2MnGa, we studied the 
evolution of the TA2 phonon anomaly with uniform 
volume-preserving tetragonal strain. 

Studying the phonon dispersions of the cubic Heusler 
structure {c/a, = 1) and its tetragonal distortion with 
c/a = 0.94, we have identified the existence of the un- 
stable mode at C = 0.34 and ( = 0.43, respectively, that 
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explain the transitions to the premartensitic and marten- 
sitic phases with decreasing temperature. The fact that 
for the two tetragonal strains the soft mode is always 
the TA2 mode but with a different wave vector suggests 
that the wavevector of the soft mode and the tetragonal 
strain are coupled. We therefore studied in more detail 
the dependence of ( upon tetragonal strain. 

The dispersion of the TA2 branch along the [110] direc- 
tion is shown in figure 8 for several values of the c/a ratio. 
An overall softening of the TA2 branch is observed with 
decreasing c/a ratio. With compressive strain (c/a < 1) 
the anomaly occurs at a larger (q than in the cubic struc- 
ture, and becomes broader and more pronounced. For 
c/a < 0.91 the entire branch is unstable. With tensile 
strain (c/a > 1) the anomaly shifts to a smaller (0 and 
becomes less pronounced. The lattice instability is com- 
pletely eliminated for c/a ^ 1.06. This is in agreement 
with the experimental observation of no superimposed 
modulations in a stress induced tetragonal phase with 
c/a = 1.18.25 

In figure 8, we have shown the TA2 branch for several 
different values of c/a. For each of these values of c/a, 
we can find the value of the wavevector (q for which the 
imaginary frequency of the TA2 branch has its minimum. 
This gives us the dependence of the modulation wavevec- 
tor upon strain and is shown in figure 8. Three regimes 
can be identified: (i) c/a < 0.91, where the entire TA2 
branch is unstable and there is no unique minimum; (ii) 
0.9 < c/a < 1.06, where the TA2 branch displays a well- 
defined dynamical instability localized at the wavevector 
^0 that depends on c/a; and (iii) c/a}^ 1.06, where the 
lattice instability is completely eliminated and the TA2 
branch is stable over the whole BZ, so that no minimum 
is defined. Only in regime (ii) can a modulated structure 
be expected to occur. 

In order to quantify these trends, in figure 10 we de- 
compose the phonon frequency u)^ into a short-range 
"normal" part u>^ and a long-range "anomalous" part ui^: 
uj^{()=ujI{C) + LoKC). The "normal" contribution, de- 
fined as co^{Q=A[l — cos{(Tr)], depends only upon short- 
range interatomic force constants and is connected to the 
local chemistry of the crystal. It corresponds to the sim- 
ple model for the transverse [110] vibrations of the [110] 
planes when they interact only with a first-neighbor in- 
terplanar force constant: K=-Am/2, where m is the mass 
associated with each plane. The long-range "anomalous" 
contribution depends on electronic screening effects 
that can be strongly q dependent. To fit the constant A 
in the definition of w^(C), we assume that for ^=1 the 
anomalous contribution is zero. 

The short-range contribution softens with decreasing 
c/a ratio, causing the overall softening observed for the 
TA2 branch. Decreasing c/a corresponds to a volume- 
preserving tetragonal distortion with a shorter c and 
larger a. Therefore smaller c/a values correspond to a 
larger [110] interplanar distance which leads to a weaker 
interplanar force constant and therefore a softer w^. 

The anomalous contribution is negative, large, and 



peaked at a critical wave vector Co, indicating a strong 
screening due to the electon-phonon coupling and Fermi 
surface nesting. Both Co and the intensity of the anomaly 
are strongly dependent upon strain. The value of Co in- 
creases with decreasing c/a, the intensity reaching its 
maximum at c/a=0.97. 

V. SUMMARY AND CONCLUSIONS 

Our ab initio study of the phonon dispersion and lat- 
tice instabilities in cubic and tetragonal Ni2MnGa can be 
summarized as follows. 

First, the parent cubic phase exhibits a marginal elastic 
instability, C « , meaning that there is almost no 
energy cost for a small uniform tetragonal distortion. 

Second, for 0.91 < c/a < 1.06 a Kohn anomaly is 
present due to (^l(x;tron-phonon coupling and Fermi sur- 
face nesting. It develops into a deep minimum corre- 
sponding to a localized dynamical instability whose wave 
vector Co is related in a one-to-one fashion to the c/a 
ratio. Therefore, tetragonal structures with c/a in this 
range are unstable towards a particular transverse (110)- 
shufHing modulation with a specific wavevector Co, the 
one-to-one relationship of the c/a ratio and the modula- 
tion period being due to the behavior of the topology of 
the Fermi surface under tetragonal distortions. 

Consequently, the energy of the cubic structure may be 
lowered by any one of a family of deformations, {e, u(Ce)}, 
consisting of a tetragonal strain, e, and a superimposed 
(llO)-shuffling modulation, u(Ce), whose periodicity, C, is 
directly related in a one-to-one fashion to the c/a ratio. 

The near-vanishing of C together with the fact that 
the energy surface as a function of c/a is remarkably flat 
permits a modulated ground state with c/a significantly 
different from 1. Indeed, the computed relation between 
c/a and the modulation wavevector is in excellent agree- 
ment with structural data on the premartensitic (c/a = 
1) and martensitic (c/a = 0.94) phases of Ni2MnGa. and 
also with the unmodulated stress-induced phase (c/a = 
1.18). 

The next step is to carry out total energy calculations 
of the complex modulated structures in order to deter- 
mine the equilibrium amplitude and phase of the lattice 
modulation and the associated energy gain. The energy 
could then be expanded in terms of symmetry invariants 
of the relevant structural degrees of freedom (strain and 
modulation) and determine the coefficients of each term 
of the expansion using ab initio calculations. In this way 
one could develop an ab-initio based finite-temperature 
statistical model of the phase transition in Ni2MnGa. 

The three phases of Ni2MnGa observed with increas- 
ing temperature can be understood as follows. At low 
temperatures, a minimum-energy modulated tetragonal 
structure is expected, consistent with the observed c/a 
= 0.94 phase. As temperature increases, the favorable 
entropy arising from the larger phase space for the low- 
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frequency phonons in the cubic structures leads to a 
phase transition to the cubic-modulated structure, with 
a modulation wavevector shifted down from that in the 
tetragonal phase in accordance with our results. Finally, 
the frequencies of the unstable phonons associated with 
the modulation are anharmonically renormalizcd at the 
highest temperatures and the observed structure is the 
cubic Heusler structure, with large fluctuating local dis- 
tortions. 

In conclusion, we have provided a first-principles mi- 
croscopic explanation for the origin of the physically 
interesting and technologically important martensitic 
phase transitions in the ferromagnetic shape-memory 
Heusler alloy Ni2MnGa. 
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TABLE I. Equilibrium lattice parameter, oo, bulk modu- 
lus, Bo, and magnetic moment, no, for Ni2MnGa in the L2i 

Hcuslor structure. Calculations using the spin-polarizod GGA 
(a-GGA), are compared with experimental data and with the 
results obtained using LSDA. 





ao(a.u.) 


Bo(Mbar) 


A»o(mb) 


0--GGA 


11.03 


155 


4.27 


LSDA 


10.74 


202 


3.92 


FLAPW-o--GGA'= 


10.98 


156 


4.09 


Exp. 


11.01" 


146'' 


4.17" 



"Ref. 27 
^Rcf. 23 
'^Ref. 18 



TABLE IV. Elastic constants, in GPa, for ferromagnetic 
Ni2MnGa in the cubic Heusler structure. Theory 1, elas- 
tic constants from stress calculation under strain. Theory 2, 
elastic constants estimated from the slope of the [CCO] acous- 
tic branches close to F: C', C44, and Cl correspond to the 
TA2, TAi, and LA modes, respectively. Room temperature 
measurements are shown for comparison. 





C'li 


('12 


G' 


('44 


Gl 


theory 1 


153±2 


148±2 


2.5±2 






theory 2 


138±9 


143±9 


-2.5±5 


100±5 


240±5 


Exp." 


152 


143 


4.5 


103 


250 


Exp.'' 


136±3 




22±2 


102±3 


222±9 



"Ref. 23 
^Ref. 24 



TABLE II. Phonon frequencies (in cm~^) computed at the 
zone-boundary point X. Results obtained using different ap- 
proximations for the exchange and correlation energy, cr-GGA 
and LSDA, are compared with the experiment. 





(T-GGA 


LSDA 


LSDA 


Exp. 


ao 


ao (cr-GGA) 


a„(LSDA) 


ao (a-GGA) 






82 


76 


70 


87 


X'l 


160 


179 


151 




X5 


174 


191 


161 






176 


195 


166 




xi 


178 


185 


159 




Xi 


191 


203 


168 






227 


249 


212 




x^ 


259 


286 


233 





TABLE III. Dependence upon the smearing parameter, a, 
of the frequency of the three Ss-modes at the critical wavevec- 
tor, qo. The ^3(1) mode is the anomalous TA2 mode, the 
other two, S3 (2) and S3 (3), are higher frequency optical 
modes. The k-point mesh needed to achieve convergence for 
each of the values of a is also reported. Frequencies are in 



<T(Ry) k-mesh S3(l) S3 (2) S3 (3) 

0.03 (6 6 6) 8 195 244 

0.01 (10 10 10) 39i 195 242 

0.005 (12 12 12) 39i 195 242 




FIG. 1. The fee L2i Heusler structure of Ni2MnGa. The 

small circles, large open circles and large filled circles repre- 
sent Ni, Mn and Ga, respectively. 
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100 i ' ' ' ' ' 

1 0.5 0.25 0.5 0.75 1 0.5 0.25 0.5 

(COO) (2CC0) (CCO) (CCO 

FIG. 2. Full phonon dispersion of ferromagnetic Ni2MnGa in the fee Heusler structure, along high symmetry lines of the 
fee BZ. Solid lines are fully converged ab initio calculations obtained using an effective electronic temperature tj^O.OlRy. The 
dashed line is the dispersion obtained using a higher effective electronic temperature, cr=0.03Ry, showing the dependence of the 
anomaly in the lowest acoustic branch along the (CCO) direction on the effective electronic temperature. The circles indicate 
the neutron scattering data from Rcfs. [ 19, 4]. The imaginary values of phonon frequencies are plotted along the negative 
frequency axis. The wave- vector coordinate C is in units of ^ . 
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0.2 0.4 0.6 0.8 1 



q=^(110) 2%/a 

FIG. 3. Partial phonon dispersion of Ni2MnGa in the fee 
Heusler structure, along the F-K-X line in the [110] diree- 
tion. The theoretical data are the same as in figure 2. The 
experimental data taken at 250 K and 270 K are shown for 
comparison. 
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X K F 




r qo K X 

[110] 

FIG. 4. Contour plots of the imaginary phonon frequencies 
in the (001) plane (a) and (iTO) plane (b) of the fee BZ. These 
contour plots show the region in q-space where the crystal is 
dynamically unstable. This occurs in a narrow valley situated 
along the [110] direction with its minimum at the critical wave 
vector qo=(0.34, 0.34, 0). 




FIG. 5. Fermi surface of minority spin bands for cubic 
Ni2MnGa. 
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300 I ' ' ■ I ' ' ' I ' ' ' I ' ' I I ' ' ' I I ' ' ' I ' ' ' I ' ' ' I ' ' I I 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

q=C(110)27t/a q=!;(0 1/a 1/c) 2n 

FIG. 6. Ab initio [110] and [Oil] phonon dispersions for fer- 
romagnetic Ni2MnGa in the tetragonal structure (c/a=0.94, 
0=11.26 a.u.) 




FIG. 7. Fermi surface of minority spin bands for tetragonal 
NiaMnGa (c/a=0.94). 
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FIG. 8. (color online only). Dispersion of the squared fre- 
quencies, for the anomalous TA2 phonon branch, computed 
for different volume- preserving tetragonal distortions, with 
c/a ranging from 0.88 to 1.06. 
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FIG. 9. Dependence upon c/a of the modulation wavevec- 
tor ^0 and the related modulation period (in unit of (110) 
interplanar distances). 
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FIG. 10. Decomposition of lo^ into "normal" (cjjjj) and 
anomalous [lo^ parts for different volume- preserving tetrag- 
onal distortions, with c/a varying from 0.91 to 1.06. Dashed 
lines are for c/a > 1, dot-dashed lines are for c/a < 1 and 
solid lines are for c/a = 1. 
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